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Abstract 

Two approaches (TW and ASvM) to derivation of integrable differential equations for 
random matrix probabilities are compared. Both methods are rewritten in such a form 
that simple and explicit relations between all TW dependent variables and r-functions of 
ASvM are found, for the example of finite size Gaussian matrices. Orthogonal function 
systems and Toda lattice are seen as the core structure of both approaches and their 
relationship. 



* e-mail: igorrumanov@math . ucdavis . edu 



1 



I. INTRODUCTION 



Random matrices were introduced in the middle of last century by Wigncr and Dyson as 
a description of complex systems that were difficult to study since they did not admit any 
exact or even approximate solutions. So the fact, thoroughly realized in the early nineties, 
is still amazing: many interesting random matrix (RM) models are closely related with and 
conveniently described by integrable systems — integrable hierarchies of partial differential 
equations (PDE). The first such result by Jimbo, Miwa, Mori and Sato^^ appeared earlier but 
remained unrecognized as an instance of a general rule for a decade. Matrix integrals, giving 
the probability distributions for eigenvalues of random matrices, proved to be r-functions 
of integrable hierarchies. The ground-breaking works of Tracy and Widom^^'^''' followed, 
which generalized and extended the results of Ref. 13 for the sine kernel to a number of 
other cases. The authors derived integrable PDE or Painleve ODE as equations satisfied by 
the probabilities for the spectrum of various random matrices as well as by their large size 
asymptotics. 

Another major approach was developed by Adler, Shiota and van Moerbeke^"^. It is 
based directly on the correspondence between matrix integrals over the spectral domains and 
T-functions — the "generating functions" of integrable PDE. Using both methods, different 
integrable PDE have been obtained not only for single matrix ensembles but also for coupled 
Hermitian matrices^'^'^'^^. In the latter case equations derived by difi:erent methods seem to 
be so dissimilar that the only indication of the existence of relations between them is that 
they are equations for the same probabilities. 

In this paper, we compare the "algebraic" approach of Adler, Shiota and van Mocrbeke 
with the "functional-theoretic" approach of Tracy and Widom and show the common structure 
behind them. This structure is (not surprisingly) the Toda lattice^^'^'^. The relation between 
these two approaches has not been heretofore completely clarified. The reason hes in the fact 
that the authors of Refs. 2, 3 used the KP hierarchy rather than one-dimensional Toda lattice 
hierarchy (1-Toda) to derive their PDE for single Hermitian matrices. This led to equations 
very difi^erent in form from those derived in Refs. 16, 17. J.Harnad^^ demonstrated for the 
special cases of Airy and Bessel kernels, how the "KP" -equation of ASvM can be obtained by 
taking the appropriate combinations of TW equations. We show that there is a more direct 
and complete correspondence of structures involving the Toda lattice hierarchy; and thus the 
procedure of Ref. 11 is essentially an implicit transition from 1-Toda to KP form of integrable 
equations related to Hermitian 1-matrix models. In this paper we restrict ourselves to the 
Hermitian one-matrix Gaussian case. Our main purpose is to reveal explicitly the underlying 
structures connecting dependent variables and equations in the different approaches. As we 
show by this example, for single matrices equations obtained by different methods (Ref. 17 vs. 
Refs. 2, 3) can essentially coincide, despite the very different ways to arise. For comparison, 
we give the implementation of the Hamad's approach in the Appendix for the case of finite 
size Gaussian matrices considered here. 

In another related work^*^, 1-Toda lattice was ingeniously used to treat the Virasoro con- 
straints of ASvM approach, which, together with an equation of Ref. 17, helped to find 
Painleve VI equation for the Jacobi ensemble. 
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II. SUMMARY OF RESULTS 

We use the ASvM approach of Hirota's bihnear identities for r-functions of the 1- dimensional 
Toda lattice hierarchy^^ and (linear) Virasoro constraints. Here the r-function has the 
meaning of complex Hermitian Gaussian matrix integral over spectral domain J"-, J = 
Ui(a2i_i, a2i) - a finite union of intervals of the real axis M with endpoints af 

where n is the random matrix size, t — {ti,t2,t3, ...) is a complex infinite vector, pt{x) — 
Q-v{x)+Y.'^^itk'-c _ f^i^Q t-deformed weight; and T„(i) will denote the corresponding integral 



over the whole W^. Our first results are the equations 

B\lnT^ = 2{2Qr,Pn-n), (1) 

Bl.Qn = (-2Bo + An)Qn - SQlPn, (2) 

B\Pr, = (2Bo + 4n)P„ - SQ^P^^ (3) 



where Qn = t^+i/Tuj Pn = ^i-\l^n^ ^-iid the operators By^ contain partial derivatives w.r.t. 
the boundary points of the intervals: 

i(i{end'points\ 

Then, analyzing the TW equations for the finite n Hermite case, we obtain simple relations 
among the above r-functions and their ratios and the auxiliary functions -u, f,t(7, appearing 
in Ref. 17. The TW approach consists in studying the kernel of a Fredholm integral operator 
K^^ whose determinant is the probability of having no eigenvalues in a subset of R, here - 
the complement J"^ of the set J: 

det(7-X„^) = ^. (4) 

where the numerator and denominator of the r.h.s. are the matrix integrals (the r functions 
rit — 0)) over J and over real axis R, respectively. The operator has kernel 

y) = Kr,{x, y)xAy), (5) 
where Xjc(x) is the characteristic function of J'^: 

if xer 
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Kn{x,y) is an integrable kernel^^, which can be expressed by Christoffel-Darboux formula for 
orthogonal functions, 

j,^^^^^^^mm^mmjf2M-)My). (6) 

^ k=0 

(f){x) — ■\Jhn-i(t>n{x), — \/bn-i(t>n-i{x), and is the coefficient of in the 

corresponding 3-term recurrence relation. 

We show that the inner product v — (0, (/ — K^)~^%1)) is related to in the following 
way: 



.^(0,(7-i^„^)-V) = ^=^. (7) 



Also, if we denote 



Qn = -u, Pn = w + vW2, (8) 

where u = (0, (/ — K:^)~^(p), w = (/ — K^)~^ip), then these functions Qn and P„ turn 
out to satisfy the same nonhnear PDE (2) and (3) as the ratios of r-functions Qn and P„, 
respectively. These facts lead to even stronger conclusions: 

Qn CnQni Pfi Qn/ 

where C„ = (2"/n!) ■ ^n/27r. So the functions u and w appearing in the resolvent kernel 
approach, are represented in terms of ratios of matrix integrals (r-functions at t = 0): 

u ^ (0, (7 - Kir4>) - a/S fl - -^^1 = ^^^^ i 1 - , (9) 



2 V n!0r J ]J (r^)^ V r^/r„ 



w 



(/-^„r V) = (^^^ - 1) ^ - 1) . (10, 



Thus, the correspondence between TW auxiliary variables and r-functions of one-dimensional 
Toda lattice hierarchy considered by Adler and van Moerbeke,^''^ becomes direct and explicit. 
The probabilities, expressed in the form of r-ratios are seen from the above formulas to be 
also the natural and convenient auxiliary coordinates. Indeed, r-functions and their ratios 
seem to provide the best description for all integrable problems. Here, in the Gaussian case, 
for instance, we have the following recurrence relations for the r-ratios: 



1 

Qn 



Pn+1^^. fill 



1 lnQ„ ^ , , , 

Qn+l = Qn { -\ + QnPn + ^), (12) 
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the first — by definition of Qn, Pn, the second follows from Eq. (1) and the first one. Thus in 
principle these ratios can be recursively determined from the simplest ones {Qq and Qi here, 
with our normalization, Qq — t( — fjS'^ dx, Qi = t//t/, = 1/2 fj2 ~^ {x — y)^dxdy). 

Remark. The 1-Toda lattice r-functions arc naturally related to the orthogonal polynomial 
structure and this is the reason for the existence of such simple relations. The orthogonal 
functions for t-deformed weight satisfy the t-dependent three-term recurrence relation 

x(j)n{x; t) = an{t)(t>n{x] t) + t) + t), 

where the coefficients depend on times and are known in terms of the Toda r-functions^'^: 

d Tn+l 
dti Tn 



r2 

' n 



These relations show that the 1-Toda time flows are in fact flows between different orthogonal 
function systems. The 1-Toda bilinear relations may also be considered as a form of orthog- 
onality relations for systems of orthogonal functions. Recall also that Christoffel-Darboux 
formula, which is the starting point of the resolvent kernel approach, is itself equivalent to 
the three-term relations for the corresponding orthogonal functions. One can see that the 
structure of orthogonal function systems is at the heart of all integrability derivations and 
connections with random matrix theory. 



III. DERIVATION OF EQUATIONS (1),(2),(3): BILINEAR IDENTITIES AND 
VIRASORO CONSTRAINTS 

We consider Hermitian random matrices of size n x n with spectral measure d/x = p{x)dx, 
p{x) — e~^^^^ such that all its moments exist. Then the probabihty that all eigenvalues are 
in J C R is 

fj„Al(x)Y\^pixi)dxi , , 

^IkA 111 py^iY^i 

where A„(a;) = ni<j(^i ~ ^i) the Vandermonde determinant. The r-function of integrable 
hierarchy appears from deformation of the above integrals: let t — (ii, ^2, ^3, ■■■) be a complex 
infinite vector of "times" and pt{x) — e~^(^)+^fc^i - the deformed weight, then the matrix 
integral 

1 /■ 

^nit)-- Al{x)l[pt{x,)dx„ (14) 

is a r-function of integrable hierarchies such as KP and one- dimensional Toda lattice hierarchy 
(see, for instance, Refs. 1 — 3). 
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Matrix integrals for coupled random matrices satisfy*^ Hirota bilinear identities for 2- 
dimensional Toda lattice (2-Toda)^^: 

/ Tn{t - [Z-'], s)Tm+l{t' + [^~'], S^^^^^^-''^^' z"-'dz = 
J z=oo 

= / TnMt,s-[z])Tr.{t',s' + [z])e^^^'^-''^>-'z'^-'--'dz (15) 

Jz=0 

The 2-Toda r-functions have two infinite sets of time parameters. There are two possible 
reductions of them to the biUnear identities for 1-matrix integral r-functions, possessing only 
one set of times (first was written in Ref. 3, second - in Ref. 6). The first one can be 
obtained by setting r{t. s) = r(t — s) and s — s' = —{t — t'), then renaming variables t — s t, 
t' — s' t' . It gives the bilinear identities^ for 1-dimensional Toda lattice hierarchy (1-Toda) 
(including the AKNS hierarchy and the one- dimensional Toda chain itself, see, e.g., Ref. 19): 

/ rn{t - [z-'])T^+^{t' + [;^-i])e^/2Er(t.-t'.>'=^n— -1^^ ^ 

J Z=CX) 

= <f r„+i(t + [;2])T^(t'-[^])e-i/^^?°(*'=-*^)^"'z"-"^-Mz (16) 

The second one is obtained by setting r(t, s) = T{t — s) and s = s', then again renaming 
t — s ^ t, t' — s' ^ t'. This gives the bilinear identities^ for the discrete KP hierarchy for 
n > m: 

J z=oo 

= (f Tn+lit + [z\)T^{t' - [z])z"-^dz = Oforn>m (17) 

The most familiar bilinear identity of the usual KP hierarchy^ corresponds to taking m — n — 1 
in (17). The integrable nonlinear PDE can be readily obtained from the above identities by 
shifting times t ^ t + a, t' ^ t — a, expanding in Taylor series in a and taking residues for 
each term of this expansion independently. 

Let us proceed in the 1-Toda setting, which, as we will see, shows most clearly the corre- 
spondence between the ASvM and the TW approaches to the derivation of nonlinear integrable 
PDE for the spectral probabilities of random matrices. 

The simplest and most important equations are obtained from (16) for m = n and m = 
n — 1 upon taking the residues of the terms linear in 2a = t — t'. For m = n one gets the 
hierarchy of equations in terms of r-functions 

k ^ ~ 

d Tn+l _ Pi{-dt)Tn ■ Pk-jjdtjTn+l ^^^^ 

dtk Tn ^ TnTn+1 

where Pk{t) are the elementary Schur polynomials defined by 
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eEr*.^^^J^p,(t),^ (19) 

fc=0 

dt = {dt^-i 1/2^^*2) l/3f?t3, th^t one can write 

oo 

rn{t T [z-^]) = 5Z;^-'pfc(T4)rn(t). (20) 

fc=0 

Equation (18) is trivial when k — 1. The simplest nontrivial equation arises from (18) when 
k — 2. It can be written as 

— ln^ = --2ln^+ — In^ +2^72^, 21a 

dt2 T„ dt{ Tn \dti Tn J dt{ 

or, introducing the new variable Q„ = r„+i/r„, as 

Now consider the 1-Toda bilinear identity for m = n — 1. The same procedure of taking 
residue for terms of order a as above now gives another infinite set of integrable PDE: 



o d'^ , ^ Pi{-dt)-Tn ■ Pk+l-i{.dt)Tn _ Pi(-4)Tn+l ■ Pk-l-i{,dt)Tn-l 

Then the simplest nontrivial equation is obtained when A; = 1, and this is the standard Toda 
equation for 1-dimensional Toda chain expressed in terms of r-functions: 

d'^ In Tn _ Tn+lTn-l , . 

dtl - Tl ■ 

Thus we obtained a closed system of three nonlinear integrable equations for three dependent 
variables, lnT„, Qn and P„: 

= QnPn, (23) 
dQn _ d^Qn ^ dHuTn ^ 

dh ~ dtj dtl ^ ^ 

dPn d^Pn d'^lnTn^ 

+ 2^^^n, (25) 



dt2 dtl dtl 

which is nothing but the well-known Toda-AKNS system (see, for instance, Ref. 15). 

Next we will use the well-known Virasoro constraints to express the time derivatives in 
the previous equations in terms of boundary derivatives, i.e. the derivatives with respect 
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to the endpoints of the intervals where the matrix eigenvalues lie. To make the paper self- 
contained, we reproduce here the derivation of Virasoro constraints for the general (5 one- 
matrix ensemble/'^'^ Although for the purposes of current paper we will need only the simplest 
case (3 = 2 corresponding to Hcrmitian ensemble, it is instructive and not more difficult to 
give the derivation for arbitrary real nonnegative j3. The Virasoro constraints arise from 
the obvious invariance of the matrix integral (the r-function) under an arbitrary change of 
integration variables. Consider infinitesimal such changes of the form Xi — > Xi-\-ex\'^^^ £ — > 0, 
A; > — 1. Then the generators Vk of such changes of coordinates make a subalgebra of Virasoro 
algebra. The change in the integrand, which will be expressed in terms of the time derivatives 
of the integral, is compensated by the corresponding change of the limits of integration, which 
in turn will be expressed in terms of the boundary derivatives of the integral. Consider the 
change in the integrand din of the matrix integral, ^'^'^ 



d{dln{x + ex^+^)) 
de 



e=0 J J 

If A; > the change of the Vandermonde determinant A„ is 



/ dinix) = T^it) ^ - j ^i{x)\{pt{x,)dx,. (14') 
J J" ^- J J" 



dA^{x + ex''+^] 
de 



also 



£=0 



de 



Using the identity 



n 

i=l 



d{d{xi + ex^i^')) 
de 



e=0 



n n k+1 _ fe+1 



i=l j>i 



X ^ X j 



\ fci+fe2=A;,fe(>0i,j=l ^ 



n n 



£=0 



^(k+i)j2^^-l[dx,. 



i=l i+1 



d 



— din ^Y.-^. din. 



i=l 



one gets (the second term on the r.h.s. is absent for k = 1) 

d{dln{x + ex^+^)) 



de 



£=0 



(26) 



E 

1=1 



fc-i 



X, 



dxi 2 dtidtk^i 



d 



1=1 



l+k 



dtu 



8 



For k = —1 the measure A^(a;) J^^'iixj does not change under translation x ^ x + e, so 
one has instead 



d{dln{x + e)) 
de 



e=Q 



^dV ^ , 9 



1=1 



din- 



For k — Q one has for the change of the measure 



d{K lii dxi{x + ex)) 



de 



e=0 



-n(n 



-l)+n] A^^{x)l[dx,, 



and one gets 

d{dln{x + ex)) 
de 



£=Q 



EdV ^ , (9 13 . ^.\,r 

i=l * 1=1 ' 



Thus one finally obtains the k — —1 constraint: 

dV 



dxi 



nt 



d 



1=2 



dti^i 



dln{x), 



the A; = constraint: 



{ i=i * 1=1 '■ ) 

and the k > constraints (again the second term on the r.h.s. is absent for k — 1): 



(27a) 



(276) 



- E ^r^^ + f E ^ + [/^^ + (1 - + + E ^^^^ ^ 



If the potential V{x) is polynomial, ^ — ^i=q9ix^ ^ then one has 



(27c) 



dxi ■^-^ Qt]j^i;^-\ 

i=l * Z=0 t+K+i 



The infinitesimal changes of integration limits are expressed in terms of boundary operators^ ^ 

d_ 
dai 



B.ri^ Y: «r^^n'- (28) 



iE{endpcnnts} 

In general, also the next important consequences of the above Virasoro constraints are needed^"^: 
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BkBiT^^ = BkViri = ViBkT'^ = ViV^ri, (29) 

the second equality holds because the operators Bk and Vi commute. 

Prom now on we consider the special case of matrices with Gaussian distribution of matrix 
elements, i.e. Gaussian potential V{x) = x^. At the locus ti = ^2 = ■■■ = 0, the Virasoro 
constraints (27a), (27b) together with identity (26) give in this case^'^: 

^Inr"^ B ^ \n.T^ 

^ 'n *-'_iJ.j.±j^ (30) 



1 



~dq~ ~ 4 + 2 ' ^^^^ 



2 



These are the only Virasoro constraints we will need here. More such relations (involving also 
the use of (27c) for A; = 1) are needed for the case of KP-hierarchy,^'^ and the derivation is 
longer. We substitute the relations (30)-(32) into (23), (24), (25) and after a simple algebra 
obtain the "boundary form" of equations (23)-(25), which is our system (l)-(3). 

If we used the KP hierarchy to derive the PDE as Adler, Shiota and Van Moerbeke did 
in Refs. 2, 3, rather than 1-Toda, we would have used just one standard integrable PDE — 
the KP equation. 



In T„ „ 9^ In Tr, . (9^ In r„ / (9^ In r, 



+ 3^:r^-4^-^ + 6^^ =0. (33) 



2 



dtf dtl dtidts V 9ti 

Then we would need more Virasoro constraints like eqs. (30)-(32) (involving also the third 
time t^, the second derivative w.r.t. t2 and the fourth derivative w.r.t. ti). Then the ASvM 
boundary KP equation would come out, for the finite n Gaussian case it reads^: 



(Bli + 8nBli + 12-B^ + 24Bo - 16B_iBi) In + 6{Bi^ In r^)" = 0. (34) 

IV. TRACY-WIDOM EQUATIONS AND r-FUNCTIONS: THE GAUSSIAN 
CASE 

We consider the same Hermitian Gaussian random matrices as above. Various large n limits 
for this case belong to important widespread universality classes: bulk limit leads to equations 
for the Sine kernel, whereas edge limit (largest eigenvalues) leads to Airy kernel. ^^'^'^ One 
studies the operator (4) with the kernel (5), (6) and = for the Hermite case. 

Also important for this theory is the resolvent kernel of K (we will use short-hand notation 
K for K;^), R{x, y), the kernel of K{I — K)~^, which is defined by the operator identity 

{I + R){I -K)^I. (35) 



10 



Then one introduces^'' the auxihary functions (denoted in Ref. 17 as Q and P, respectively) 

q{x- J)^{I- i^)-V(^), P{x; J)^{I- K)-'^{x), (36) 
and auxihary inner products, which are functions of only the endpoints Ofc of J, 

u = {q, 0Xjc), V = {q, ipXJ-) = (P, 0XjO, f« = (p, '4^Xj-)- (37) 
Then the resolvent kernel R{x, y) is^'^'^'^ 

rji N g(x; J)p(|/; J) - p{x- J)q{y- J) 

R{x,y)^ , {x,y ^ J , X i^y) 38) 

x-y 

and R{x,x) = p{x; J)q'{x; J) — q{x; J)p'{x; J). The TW equations for 1-matrix case are the 
equations for the functions of the endpoints — R{aj, a^), 

qk^q{ak;J), pk^p{ak]J), 

and the auxiliary functions m, t;, w. Among them there are universal equations valid for any 
Hermitian 1-matrix unitarily invariant model. They read^^ 

|^ = (-l)^i?(a„a,)gfc, (39) 
oak 

l)^R{aj,ak)pk, (40) 



% 



for j ^ k and 



i?,,^i?(a„a,) = ^^^^^, (41) 

Qj — Ofe 



Rk = R{ak,ak) ^ Pk^ Qk^, (42) 

i-^M, (43) 



dw 



= (-!)'», (44) 

The other equations are not universal, their particular form depends on the potential V{x). 
For the finite n Gaussian (Hermite) case they are^^ 

1^ = -ajqj + [yf^ - 2u)pj - ^(-l)''i?(a,-, ak)qk, (46) 
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^ = ajPj - ( + 2w)qj - ^{-\fR{aj, ak)pk, (47) 

along with 

R{aj,aj) = -2ajqjPj + {V2^-2u)p'^j + {V2^ + 2w)q] + ^{-lfR{aj,ak){qjPk-Pjqk): (48) 

= -2g,p, - 5:(-l)'=i?(a„ a.)^ (49) 



It follows from (39)-(41), (46) and (47) that 

d 

and the last equation together with (43), (45) gives the first integral obtained in Ref. 17: 

2 '^{-lyqjpj = -{2u - V2n){2w + V2n) - 2n, 
j 

or, if we denote V2n — 2u — 2(5„, 2w + V2n — 2P„, as in formula (8), we get 

J](-l)%p, = 2g„P„-n, (50) 
j 

Prom (39)-(41), (46) and (47) it follows also that 

+ - 2ii)p,-, (51) 

k ^ 

Then consider the expression (using (43), (51), (8) and (50)) 

= 2 |^afc(-l)'=5^_2g^^(_l)fe^^pj ^ _2 ^ - 4g„(2g„P„ - n), 

which is nothing but 
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B\Qn = 2(-Bo + 2n)Qn - SQlPn, (2') 
in terms of the earlier introduced Qn- Completely analogously one can obtain 

B\Pn = 2(Bo + 2n)Pn - SQnPl (3') 

in terms of P„ introduced in (8). 

Now let us prove the last equation 

B\\nT^ = 2{2Qr,Pn-n) (!') 

with the same Qn and P„ defined by (8). Recall formula (4) for the Predholm determinant as 
a ratio of r-functions. On the one hand, it follows from (4) that 

^-1 l^^n = E ^ l^^n = E ^ l^det(/ - K^), (53) 



on the other hand, 



^ dak k 



k k 



= Tr{I - Ki)-^ Yl ^Kn{x, y)xAy; a) = - ^(-1)'^^- (54) 

Thus, we have 

S_ilnr„^ = -5](-l)^'i?,. (55) 

k 

Taking expression (48) for R). = R{ak, Ofe) we get 

Y,{-^fRk = -2 Y.^-l)\auqm - QnPl - Pnql) - E {-iy^''Rjk{p,qk - QjPk), 

k k jjkytj 

the last double sum above, however, is equal to zero. By (43), (45), the last expression is 

J](-l)^'i?, = -2 J](-l)^a,g,p, + E 1^ - E (^6) 

k k k k 

Another auxihary formula we need follows from (51), (52), (43) and (45): 



E(-l)'«^^^ = '^Y.^-lfak{QnPl-Pnql) = 2E«M Qn'^ + Pnf^ = 2B,{Q^P^ 
k,j k k \ "'^ "'^ J 

(57) 
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Using (56), (57) and (50), we derive from (55): 



fc j k 

and, using (2'), (3'), this is exactly 

B!ilnT„^ = 2(2g„P„-n). 

Equation (1') with TW variables is proven. Thus, we have shown that Toda lattice stands 
behind TW derivations. 



V. TW VARIABLES IN TERMS OF r-FUNCTIONS 



The TW u and w variables for Gaussian case turn out to be directly related to the ratios of 
Toda lattice r-functions Tn+i/rn and Tn-i/rn- Moreover, if we compare formulas (1'), (44) 
and (50), we see that in the Gaussian case, the equation holds: 

B\ lnr„^ = 2(2g„P„ - n) = 2B.iv, (58) 

or 

B-i(B-i\nT^ -2v) = 0. (59) 

Besides, the equations (58) and (55), (56) mean, respectively, also the following relations 
among r-functions and TW variables in the Gaussian case: 

QnPn = QnPn, (60) 

B-i In = 2Bov + 2PnB-iQn - 2g„B_iP„. (61) 
Then from (2), (3) and (60) it follows immediately that 



Bl,Q„ + 2B,Qn Bl,Q„, + 2BoQn Bl,P„ -2B^Pn Bl,P„,-2BoPn ^ , 
7^ = ^ = p = 5 = 4(n - 2g„P„). 

Vn r^n 

(62) 

From equations (62) one can readily obtain by simple algebra: 

B^liPnB^lQn - QnB-lPn) = -2B^{QnPn) (63) 

as well as the identical equation for the corresponding "tilde" -quantities. Since the right-hand 
side of equation (63) is equal to the one for "tilde" -quantities, the corresponding left-hand 
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sides are equal too. Now we can apply the same argument as was used in Ref. 17 to obtain the 
crucial first integral (50) in this finite n Hermite (Gaussian) case, this time for the function 

Fn = PnB-iQn - QnB-lPn = PnB-lQn - QuB-lPn- (64) 

The argument consists in that the last equality in (64) is true because clearly both sides tend 
to zero as all their arguments tend to oo, and by (63) and (62) both are translationally 
invariant w.r.t. any vector of the form (a, a, a). Then by (62) we can rewrite (64) as 

B_i In Qn - B_i In Pn = B_i In g„ - B_i In P„ (65) 
Besides, by (62) we also have 

B_i In g„ + B_i In P„ = B_i In g„ + B^i In P„, 
and so the simple relations emerge, 

-B_ilnQ„ = -B_ilnQ„, (66) 



-B_ilnP„ = -B_ilng„. (67) 

Now by applying the above argument once again, the announced result Qn = CnQn, Pn — 
Qn/Cn with some constant C„ follows. The constant C„ is possible to determine from the 
recursion, by definition of r-ratios Qn, Pn- 

1 



Qn-l 

and the limit when J becomes the whole real line. In this limit both u and w tend to zero, 
and the limit of QnPn is just n/2. Then also Qn/Qn-i = n/2 (this is as it should be because 
QnPn = ^^^^ limit) and, comparing with the well-known value of the size n GUE 

matrix integral over the whole space (see, e.g., Ref. 14), 

n/2 n-1 

^n(0) = ^;^^ll^\, (68) 
j=i 



we find that 



On the other hand, equation (58) is also possible to integrate once by the same argument 
that in the limit of all Oj — > oo both functions v and P_i lnT„ tend to zero. So by translation 
invariance they must be equal: 

v = B.,lnT^/2. (7) 
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VI. CONCLUSION 



A clear connection is exposed between the structures of equations derived by different — 
Tracy- Widom and Adler-Shiota-van Moerbeke — methods. The structure of universal equa- 
tions arising from resolvent kernel analysis is identified with that of 1-Toda-AKNS integrable 
hierarchy. We have shown how, in the Gaussian case, to express the auxihary coordinates, 
related to the resolvent kernel in the first approach, in terms of ratios of r-functions, the 
natural coordinates of the second approach. The last ratios have here also a simple meaning 
of probabilities. The close relation between Toda lattice hierarchies and systems of orthogonal 
functions is crucial for this connection. 

Although our explicit results have been obtained for single Gaussian matrix case, every- 
thing we have done here can as well be implemented for other one-matrix Hermitian unitarily 
invariant ensembles. Also it is well-known that the (more universal) equations for various 
large n limits can be readily obtained from the ones for finite n (see, e.g., Ref. 17 or Ref. 
3). Moreover, our results pave the way to finding some similar connections for the case of 
coupled Hermitian matrices. Such relations then may facilitate the more difficult analysis of 
probabilities for the spectrum of coupled random matrices, with many possible applications 
to random growth processes and other nonstationary nonlinear dynamics phenomena. This 
work is in progress. 
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VII. APPENDIX: J.HARNAD'S PROCEDURE - FINITE n GAUSSIAN CASE 

Following notations of Harnad,^^ introduce new TW variables: 

X2j = 2iq2j, X2j+i = 2q2j+i, {Al) 

Y2j = ip2j, y2j+i = P2j+1- {A2) 

Then, in terms of these variables, the first derivatives of r-functions are the Poisson commut- 
ing Hamiltonians Gf for a finite integrable dynamical system^^ with canonical coordinates 
{Xj, Yj, u, w} in our case. 
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Gf=^-^ = {-iy^'R{a,,a,) 



the relevant TW equations are 

— — —^1^ ^k, 

oak CLj - a-k 

_ .in^jYk-YjXk^ 



for j ^ k and 



1/2^^ ^F,, 



f) Y \^ V V Y 

^ = -a,-X,- + 2(v^ - 2i.)y,- + 1/2 y ^^ii^-^l^Xfe, 



2m 

1^ = a,Y, - 1/2{V2^ + 2w)X, + 1/2 V SklM^y,. 



Let 



2m 



2m 2m 2 m 

.?■=! j=l j=l 

Then the first integral (50) is 



2V2n(M -w)+ Auw - Sq = 0. 
We express by the above formulas (let u = \/2n — 2u, w = \/2n + 2w): 

2m 

B.i Inr^ = R^ = Y^Gf ^ -Si + uPo + l/4wQo, 
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2m 



Bo In r;f = i?f = ^ a^Gf = -S^ + uPi + l/AwQi - 1/8 - r.Xfc)^ = 

= -^2 + uP^ + l/4wgi - l/AQoPo + l/4S'o^ (^8) 

we need also 

2m 

i3i In = i?f = ^ afCf = -^3 + UP2 + l/4wg2 - 1/8 ^^(a,- + a^) (X.F, - Y^Xj^f = 

= -,^3 + m + l/4«;g2 - l/KQoPi + Qi^'o) + l/'iSoSu (A9) 

because + - I^'^a^)' = 2 T.j,k ^'jiXfY,^ + Yf^l - 2X^XkY^Yf,) = 2(QiPo + 

Qo-Pi ~2S'oS'i). We now find from (A7)-(A9) all the terms appearing in the ASvM boundary- 
KP equation (34). Applying the operators Bk to (A7)-(A9), we use the Poisson commutativity^^ 
of {G'j^}j=i,...,2TO- This property means that the derivatives of Rf with respect to the a^-'s are 
just given by their explicit dependence on these parameters. ^-"^ This greatly simplifies the 
necessary differentiations so we obtain: 

B\ lnr„^ = B_^R^ = -So, {AlO) 
B_iBi Inr^ = B_iR^ = -3S2 + 2mPi + l/2wQi - l/2QoPo + l/^S^ {All) 
Bl In ri = BqR^ = -2S2 + uPi + l/iwQi, {A12) 

B\\nri=B\{B\\nri)=B\{B.,R^) = -B\So = "E^E^E^.^.' = 

i I j 

= - E[2^ • 21^- (a,F, - wX,/2) - w/2 ■ 2X,{-a,X, + 2uY,)] - E[41^' • E ^'/^ + E 

j j i i 

= E[-2 E ^i^' - 4%^^' - + ^^i^^XjYj] = 

= -2Q0P0 - AuPi - wQi + AuwSo. {A13) 
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Now we substitute the expressions (A8), (A10)-(A13) into the left-hand side of equation (34) 
and obtain 

{B\ + SnB\ + l2Bl + 24Bo - 16-B_iBi) Inr^ + 6(Sii lnT„)2 = 
= [-2QqPq - AuPi - wQi + iuwSo] - 8nSo + 12[-252 + uPi + l/4wQi]+ 



2A[-S2+uPi + l/'^wQi-l/4QoPo+l/4S^]-16[-3S2+2uP^ + l/2wQi-l/2QoPo+l/2S^]+6S^ 



= 4:So[{V2^ - 2u)(V2n + 2w) - 2n + Sq] = 4:So[-2V2n{u - w) - 4uw + Sq] = 0, 

by the TW first integral (A6). Thus, the ASvM boundary-KP equation for the finite n 
Gaussian case (34) is proved by expressing the derivatives of the r;^-function in terms of the 
TW variables, in the way of Ref. 11. 
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